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Abstract 

By extending the approach of Berry and Robnik, the limiting level spacing 
distribution of a system consisting of infinitely many independent components 
is investigated. The limiting level spacing distribution is characterized by a 
single monotonically increasing function p-(S) of the level spacing S. Three 
cases are distinguished: (i) Poissonian if /2(+oo) = 0, (ii) Poissonian for large 
S, but possibly not for small S if < /2(+oo) < 1, and (iii) sub-Poissonian if 
/2(+oo) = 1. This implies that, even when energydevel distributions of indi- 
vidual components are statistically independent, non-Poissonian level spacing 

distributions are possible. 

PACS number (s): 05.45.Mt, 03.65.Sq 
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I. INTRODUCTION 



For a bounded quantum system in the semiclassical limit, statistical properties of energy 
levels at a given energy have been intensively studied. Universal behaviors are found in 
statistics of unfolded energy levels [1-3], which are sequence of numbers uniquely determined 
by the energy levels using the mean level density obtained from the Thomas- Fermi rule (see 
Appendix A). For quantum systems whose classical counterparts are strongly chaotic, the 
unfolded level statistics are well characterized by the random matrix theory which gives level 
spacing distribution obeying the Wigner distribution [4,5]. 

For quantum systems whose classical counterparts are integrable(such systems will be 
referred to as classically integrable systems), level statistics were theoretically studied by 
Berry- Tabor [2] , Sinai [6] , Molchanov [7] , Bleher [8] , Connors and Keating [9] , and have been 
the subject of many numerical investigations. Still its mechanism is not well understood, 
the level spacing distributions are believed to be the Poisson (exponential) distribution for 
generic classically integrable systems. 

As suggested e.g., by Hannay (see the discussion of [10]), one possible explanation would 
be as follows: For an integrable system of / degrees-of-freedom, almost every orbit is generi- 
cally confined in each inherent torus, and whole region in the phase space is densely covered 
by invariant tori as suggested by the Liouville- Arnold theorem [11]. In other words, the 
phase space of the integrable system consists of infinitely many tori which have ] infinites- 
imal volumes in Liouville measure. Then, the energy level sequence of the whole system 
is a superposition of sub-sequences which are contributed from those regions. Therefore, if 
the mean level spacing of each independent subset is large, one would expect the Poisson 
distribution as a result of the law of small numbers [12]. 

As shown by the pioneering work of Berry and Robnik [10], this would be the case for the 
nearly integrable systems consisting of large number of independent components. Quantum 
energy eigenstates are considered to be 'a superposition' of various classical trajectories, 
which are connected by the tunneling effect. In the semiclassical limit, because of the 
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suppression of the tunneling, energy eigenstates are expected to be localized on individual 
region in phase space explored by a typical trajectory of classical system [14], such as tori 
and chaotic regions, and to form independent components. The level spacing distribution is, 
then, regarded as a product of superposition of energy levels contributed respectively from 
those independent components. Based on this view, Berry and Robnik observed that the 
level spacing distribution of the system consisting of a single regular component, described 
by a Poisson distribution, and N chaotic components, described by the Wigner distributions 
of equal strength, approaches the Poisson distribution when N — > +00. 

The localization of energy eigenfunctions onto individual phase space structures was 
studied by Berry as a semiclassical behavior of the Wigner function [14]. For the classically 
integrable system, it is shown that the Wigner function tends to delta function on tori 
in the semi-classical limit [15]. Such localization phenomenon of eigenfunctions has been 
observed numerically for several systems [16-18]. For example in Ref. [18], it was shown 
that the Husimi distribution function of each energy eigenstate is well localized either on the 
classical chaotic regions or on the regular regions filled with classical tori and that the energy 
eigenvalues corresponding to eigenfunctions supported by distinct phase-space regions obey 
distinct statistics. 

Those works seem to imply that the existence of infinitely many independent components 
is an essential ingredient of the appearance of Poissonian level spacing distribution and that 
the semiclassical limit is one of the mechanisms providing infinitely many components. How- 
ever, in some classically integrable/nearly-integrable systems [2,8,20,21], where one might 
expect infinitely many components, deviations from the Poissonian distribution have been 
observed. Therefore, it is interesting to explore the consequences only of the existence of 
infinitely many independent components. 

In this paper, along the line of thoughts of Berry and Robnik [10], we investigate the dis- 
tribution of unfolded energy levels when energy levels consist of infinitely many independent 
components, and show the possibility of deviation from the Poisson distribution. Hereafter, 
for the sake of simplicity, unfolded energy levels will be referred to as energy levels or levels. 
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We briefly review the Berry- Robnik theory [10]. It relates the statistics of the energy 
level distribution to the phase-space geometry by assuming that the sequence of the energy 
spectrum is given by the superposition of statistically independent subspectra, which are 
contributed respectively from eigenfunctions localized onto the invariant regions in phase 
space. Formation of such independent subspectra is a consequence of the condensation 
of energy eigenfunctions on disjoint regions in the classical phase space and of the lack 
of mutual overlap between their eigenfunctions, and, thus, can only be expected in the 
semi-classical limit where the Planck constant tends to zero, h — > 0. This mechanism is 
sometimes referred to as the principle of uniform semi-classical condensation of eigenstates 
[13,26], which is based on an implicit state by Berry [14]. The validity of this assumption 
is checked numerically in the semi-classical (high energy) region where the Planck volume is 
much smaller than the phase volume of each invariant region [16,17,22]. 

In the Berry- Robnik approach [10], the overall level spacing distribution is derived along 
a line of mathematical framework by Mehta [4], as follows: Consider a system whose classical 
phase space is decomposed into N disjoint regions. The Liouville measures of these regions 
are denoted by p, L {i = 1, 2, 3, • • • , N) which satisfy Yh=\Pi — 1- Let E(S) be the gap 
distribution function which stands for the probability that an interval (0, S) contains no 
level. E(S) is expressed by the level spacing distribution P(S) as follows; 



When the entire sequence of energy levels is a product of statistically independent superpo- 
sition of N sub- sequences, E(S; N) is decomposed into those of sub- sequences, Ei(S; pi), 



In terms of the normalized level spacing distribution Pi(S;pi) of sub-sequence, Ei(S;pi) is 
given by 







i=i 




(1.3) 
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and Pi(S] is assumed to satisfy [10] 

l°° S-pi{S- Pi )dS = -. (1.4) 
Jo Pi 

Note that the spectral components are not always unfolded automatically in general even 

when the total spectrum is unfolded. However, in the sufficient small interval [e, e + Ae] , 

each spectral component obeys a same scaling law (see Appendix A) and thus is unfolded 

automatically by an overall unfolding procedure. Equations (1.2) and (1.4) relate the level 

statistics in the semiclassical limit with the phase-space geometry. 

In most general cases, the level spacing distribution might be singular. In such a case, 

it is convenient to use its cumulative distribution function pf, 

Pi(S) = / pi(x;pi)dx. (1.5) 
Jo 

The corresponding quantity of the overall level spacing distribution is 

M(S;N) = [ P(x;N)dx, (1.6) 
Jo 

where P(S; N) is the level spacing distribution function corresponding to E(S; N). 
In addition to equations (1.2) and (1.4), we assume two conditions for the statistical 
weights: 

• Assumption (i): The statistical weights of independent regions uniformly vanishes in 
the limit of infinitely many regions; 

maxpj^O as N — > +oo. (1.7) 

i 

• Assumption (ii): The weighted mean of the cumulative distribution of energy spacing, 

N 

p( P ;N) = J2pMp), (1-8) 



i=i 



converges in iV — > +oo to p(p) 



hm p(p;N)=p(p). (1.9) 

N—>+oo 



The limit is uniform on each closed interval: < p < S. 



In the Berry-Robnik theory, the statistical weights of individual components are related to 
the phase volumes of the corresponding invariant regions. This relation is satisfactory if 
the Thomas-Fermi rule for the phase space fractions still holds in general. Here we do not 
specify their physical meaning and deal with them as parameters. 

Under the assumptions (i) and (ii), eqs.(1.2) and (1.4) lead to the overall level spacing 
distribution whose cumulative distribution function is given by the following formula in the 
limit of N — > +00, 



where the convergence is in the sense of the weak limit. When the level spacing distributions 
of individual components are sparse enough, one may expect Ji = and the level spacing 
distribution of the whole energy sequence reduces to the Poisson distribution, 



In general, one may expect /2 7^ which corresponds to a certain accumulation of the levels 
of individual components. 

In the following sections, the above statement is proved and the limiting level spacing 
distributions are classified into three classes. One of them is the Poisson distribution as 
discussed in the original work by Berry and Robnik [10]. The others are not Poissonian. 
We give examples of non-Poissonian limiting level spacing distributions in section III. In the 
concluding section, we discuss some relations between our results and other related works. 



In this section, starting from eqs.(1.2) and (1.4), and the assumptions (i) and (ii) in- 
troduced in the previous section, we show that, in the limit of infinitely many components 



M p (S) = l-(l-ji(S))exp 



Jo 



(1.10) 





II. LIMITING LEVEL SPACING DISTRIBUTION 



A. Derivation of the limiting level spacing distribution 
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iV — > +00, the level spacing distribution converges weakly to the distribution with the 
cumulative distribution function: 



M p (S) = l-(l-fi(S))exp 



fS 

— (1 — p(a)) da 
Jo 



(2.1) 



According to Helly's theorem [27,12], this is equivalent to show that the cumulative distri- 
bution function M(S; N) converges to M^S). The convergence is shown as follows. 

Following the procedure by Mehta(see appendix A. 2 of Ref. [4]), we rewrite the gap 
distribution function E(S; N) in terms of the cumulative level spacing distribution functions 
Pi(S) of independent components: 

N r ,4-™ n r r S 



E(S;N) = H pi / da(l-pi(a)) = ]J 1 - Pi da (1 - p t (a)) 
i=i L Js - 1 i=i L Jo 



(2.2) 



The second equality follows from Eq.(1.4), integration by parts and linio-^+oo a (1 — Pi(a)) = 
0, which results from the existence of the average. Then, the overall cumulative level spacing 
distribution function M(S; N) is given by 



M(S; N) = 1 + -^E(S; N) 
= 1-E(S;N)J2 



Pi - PiPi(S) 



(2.3) 



7[ 1 ~Pi So da (I - Pi(a)) 
First we consider the behavior of E(S;N). Since the convergence of Y,f=i PiPi(cr) — > p,(a) 
for N — > +00 is uniform on each interval a G [0, S] by Assumption (ii), 



A? 



log E(S;N) = $>g 



i=l 



fS 

1 - pi da {I - Pi{a)) 
Jo 



N 



= -E 



i=l L 
r5 



Pl / rf ( r(l-^( r))+O(p2) 



- / da[l-p{a-N)] + Y.O{p>) 

Jo T 
— > — / dcr [1 — p(a)] as iV — > +00, 

JO 



(2.4) 
(2.5) 



where we have used |/Xj(cr)| < 1, log(l + e) = e + 0(e 2 ) in e <C 1, and the following property 
obtained from Assumption (i) 



N N 

I 0(p^)\ < C • max pj ■ Pi = C ■ max pj — > as iV — > +00. (2.6) 
i=i i=i 

with C a positive constant. The N — > +00 limit of the sum in the right-hand side of equation 
(2.3) can be calculated in a similar way. Indeed, as 1/(1 — e) = 1 + 0(e) in e <C 1, one has 



Pi - PiPi(S) 

Ul- Pi ^da{l-pn{a)) 



N 



N 



i-Ep^) + Eo(p, 2 ; 



i=\ 



i=l 



1 - p(S) as iV -> +00. 



(2.7) 
(2.8) 



Therefore, we have the desired result: 



lim M(S; N) = MJS) = 1 - (1 - p(S)) exp 

N— »oo 



- I" (l-fl(a))da 
Jo 



(2.9) 



B. Properties of the limiting level spacing distribution 

Since Pi(S) is monotonically increasing and < pi(S) < 1, p(S) has the same properties. 
Then, 1 — p(S) > for any S > and one has 
1 f s 

- da(l-Ji(a)) — >l-p(+oo) as S -> +00. (2.10) 
o Jo 

The limit classifies the following three cases: 

• Case 1, p(+oo) = 0: The limiting level spacing distribution is the Poisson distribution. 
Note that this condition is equivalent to p(S) = for y S because p(S) is monotonically 
increasing. 

• Case 2, < p(+oo) < 1: For large value of S, the limiting level spacing distribution 
is well approximated by the Poisson distribution, while, for small value of S, it may 
deviates from the Poisson distribution. 

• Case 3, p(+oo) = 1: The limiting level spacing distribution deviates from the Poisson 
distribution for V S, in such a way that the cumulative distribution function approaches 
1 as S — > +00 more slowly than does the Poisson distribution. This case will be referred 
to as a sub-Poisson distribution. 



One has Case 1 if the individual level spacing distributions are derived from scaled distri- 
bution functions fa as 

rS 



Pi(S) = pi I fi (pix) dx, 
Jo 



(2.11) 



where fa satisfy 



r+oo r+oo 

/ fi(x)dx = 1, / xfi(x)dx = 1 
Jo Jo 

and are uniformly bounded by a positive constant D: \fa(S)\ < D (1 < i < N and S > 0) 
Indeed, one then has 

\KS;N)\<J2p-[ \fi(Pix)\dx 

N N 

< DSY.pl < DSmaxpiYPi^Q^rtS)- 



(2.12) 



This includes the case studied by Berry and Robnik [10], where the level spacing distribution 
is a superposition of a single regular component and N equivalent chaotic components, and 
the latter is expressed by the product of the scaled distributions as in Eq.(1.2). Indeed, one 
has 

N 

E BR (S; N) = exp (-p S) ]J E~ R (S; Pi ) , (2.13) 

i=i 

where the statistical weights are p« = and the individual level spacing distributions fa 
corresponding to the gap distributions E™ lGNER (S; pi) are given by 



fa(x) 



irx 



exp 



7T . 



(2.14) 



In addition, this would be the case when the system consists of identical iV components 
where the level spacing distribution is described by a scaled form as Eq.(2.11). Such a case 
is expected when there is a symmetry such as the regular polygonal billiards. 

We remark that, when the limiting function p(S) is differentiate, the asymptotic level 
spacing distribution admits the following density: 



P P (S)= (1 - p(S)Y + fl'(S) 



exp 



Jo 



/x(cr)) da 



(2-15) 
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III. EXAMPLE 



As an example of the deviation from the Poisson distribution, we study the quantum 
systems whose energy levels are described by using positive integer numbers, m and i as 
follows, 

£m,i = m 2 + a- i 2 , (3.1) 

where a is the system parameter. Such energy levels are given, for instance, by the rectan- 
gular billiard system whose aspect ratio of two sides is characterized by a [19,20]. In this 
paper, the unfolding transformation of the energy levels {e mi j} — > {e mt i} (see Appendix A) 
is done by using the leading Weyl term of the cumulative mean number of energy levels, 

7T 

€m,i = ^(^m,i) ~j l = ^m,i- (3-2) 

For a given energy interval [e, e + Ae] , i or m can be regarded as an index which classifies 
energy levels into spectral components. In this paper, i — 1, 2, 3, • • • , N, 



N 



4^(1 + 7)6 - 7T 



air 



(3.3) 



is adopted for classification where 7 = Ae/e, and [x] stands for the maximum integer which 
does not exceed x. The relative weight of each component, pi(i — 1, 2, 3, • • • , JV), is given by 



4(1+7) , ft I 1-ai 2 ' I r> ( 1 \ if 
^(1+1) V 1 + + U \W) 11 



As easily seen, p« satisfies the assumption (i); 



< i < N. 

1+7 — — 



(3.4) 



4 / r 1 / 1 \ _i / 2 / 1 



IT 



Note that the limit of infinitely many components, iV — > +00, corresponds to the high energy 
limit, e — > +00 (see Eq.(3.3)), which is equivalent to the semiclassical limit. In this limit, 
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the statistical weight of each sub-spectrum becomes sparse, since each element of p(S;N), 
PiPi(S), tends to zero: piPi(S) < maxj pj — > 0. 

In the billiard system, each spectral component obeys a same scaling law (see Appendix 
A), and thus is unfolded automatically by an overall unfolding procedure. 

Figures l(a)-l(c) show numerical results of the level spacing distribution P(S) for three 
values of a. In case that a is far from rational, P(S) is well approximated by the Poisson 
distribution (Fig. 1(a)). While in case that a is close to a rational, P(S) shows large deviation 
from the Poisson distribution (Figs. 1(b) and 1(c)). When q is a rational expressed as 
a = p/q, where p and q are coprime positive integers, P(S) is not smooth and becomes a 
sum of delta functions in the limit of e — > +00 [2,9,20], which are separated by a same step 
A, 

A 7 



4^/pq 

Figures 2(a)-2(c) show the function | log (1 — M(S; N))\ for three values of a corre- 
sponding to figures l(a)-l(c), respectively. The dotted line in each figure corresponds to the 
Poisson distribution, M(S) = 1 — exp(-S). As shown in figures 2(b) and 2(c), the distribu- 
tion function for small value of S clearly deviates from the Poisson distribution. However, 
for large value of S, it approaches a line whose slope is l(see the dashed line), and thus, the 
level spacings for large value of S obey the Poisson distribution. 

In order to compare the non-Poissonian distributions and the classification given in the 
previous section, we consider, 



D(S;N) = ~]og 



1 - [ S (1- M(a;N))da 
Jo 



(3.6) 



When N — > +00, D(S; N) approaches ^ j s (l — p(a))da and this function distinguishes the 
three cases as follows: In Case 1 i.e., where the level spacing obeys the Poisson distribution, 
liniAr^+oo D(S; N) — 1. In Case2, liniAr^+oo D(S; N) — > c as S — > +00 (c^ 1), and in Case 
3, where the sub-Poisson distribution is expected, limAf^ +00 D(S; N) — > as S — > +00. 

Figure 3 shows D(S; N) for different values of N. From this, one can think that D(S; N) 
for N = 61905, S < 10 well approximates lim^+oo D(S; N). 
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Figure 4 shows D(S;N) for the three values of a corresponding to figures l(a)-l(c), 
respectively. In case that P(S) is well characterized by the Poisson distribution (figure 
1(a)), the corresponding function D(S;N) agrees with 1. In case that P(S) deviates from 
the Poisson distribution ( figures 1(b) and 1(c) ), D(S; N) approaches a number less than 1 
for S — > +oo. Therefore these results correspond to the Case 2. In this model, we have not 
yet observed the clear evidence of Case 3. Such a case is expected when there is stronger 
accumulation of the energy levels of individual components. 

IV. CONCLUSION AND DISCUSSION 

In this paper, along the line of thoughts of Berry and Robnik [10] , we investigated the level 
spacing distribution of systems with infinitely many independent components and discussed 
its deviations from the Poisson distribution. In the semiclassical limit, reflecting infinitely 
fine classical phase space structure, individual energy eigenfunctions are expected to be well 
localized in the phase space and give independent contribution to the level statistics. Keeping 
this expectation in mind, we considered a situation where the system consists of infinitely 
many components and each of them gives an infinitesimal contribution. And by applying 
the arguments of Mehta, and Berry and Robnik, the limiting level spacing distribution was 
obtained which is described by a single monotonically increasing function p,(S) of the level 
spacing S. The limiting distribution is classified into three cases; Case 1: Poissonian if 
/2(+oo) = 0, Case 2: Poissonian for large S, but possibly not for small S if < p,(+oo) < 1, 
and Case 3: sub- Poissonian if /x(+oo) = 1. Thus, even when energy levels of individual 
components are statistically independent, non-Poisson level spacing distribution is possible. 

Note that the singular level spacing distribution can be taken into account in terms of 
non-smooth cumulative distribution M^. Such a singularity is expected when there is strong 
accumulation of the energy levels of individual components. For certain class of systems, 
such accumulation is observable. One example is shown in section III where the results 
show clear evidence of Case 2. Another example is studied by Shnirelman [21], Chirikov 
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and Shepelyansky [23], and Frahm and Shepelyansky [24] for a certain type of systems 
which contain quasi- degeneracy result from inherent symmetry (time reversibility). As is 
well known, the existence of quasi- degeneracy leads to the sharp Shnirelman peak at small 
level spacings. 

One of the interesting features of the level statistics is the level clustering which is 
described by a nonvanishing value of the level spacing distribution function at S — 0. Level 
clustering is expected for integrable systems or mixed systems, but not for strongly chaotic 
systems due to the level repulsion. For certain class of systems, rigorous results are available; 
Molchanov [7] analyzed the energy levels of a one-dimensional Schrodinger operator with 
random potential and showed that the level clustering arose from the localization of the 
eigenfunction in the semiclassical limit. Minami also analyzed a one-dimensional Schrodinger 
operator with 5-potentials and reported the same result [28]. The limiting level spacing 
distribution obtained in this paper possesses level clustering property. Indeed, when p>(S) is 
differentiate, the level spacing distribution function (2.15) has nonvanishing value at S — 0: 



It is also interesting to extend the Berry-Robnik distribution (2.13) for the level statistics 
of the nearly integrable system with two degree-of-freedom. The classical phase space of this 
system consists of regular and chaotic regions. Since the regular regions corresponding to 
the system consist of infinitely many independent regions, original proposal for the gap 
distribution by Berry and Robnik would be replaced by 



where Ef MT (S;pi) is the gap distribution function of the ith chaotic component obtained 
from the random matrix theory. The above distribution formula is classified into cases; 
Case 1', p(+oo) = 0: Berry-Robnik distribution, Case 2', < /2(+oo) < 1: Berry-Robnik 
distribution for large S, but possibly not for small S, and Case 3', p,(+oo) = 1: A distribution 
function obtained by the superposition of spectral components obeying the sub-Poisson 
statistics and the Random matrix theory. From this classification, one can see that the new 



P A (0) = l+/i'(0)>0. 




(4.1) 
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formula (4.1) admits deviations from the Berry- Robnik distribution when /x(+oo) 7^ 0. 

For nearly integrable systems with two degree-of-freedom, Prosen and Robnik showed 
numerically that the Berry-Robnik formula (2.13) well approximates the level spacing dis- 
tributions in the high energy region that is called the Berry-Robnik regime [29], while it 
deviates in the low energy region and approximates the Brody distribution. They studied 
this behavior in terms of a fractional power dependence of the spacing distribution near the 
origin at S — 0, which could be attributed to the localization properties of eigenstates on 
chaotic components [30,31]. From the above classification, one can see that the condition 
p(+oo) = should be satisfied in the Berry-Robnik regime. While Case 2' and Case 3' in 
the above classification might propose another possibilities. When the spectral components 
corresponding to regular regions show strong accumulation, the level spacing statistics obeys 
the distribution formula (4.1) with < p,(+oo) < 1, and shows deviations from the Berry- 
Robnik distribution. Such possibilities will be studied elsewhere. 
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APPENDIX A: UNFOLDING OF SPECTRUM 

The unfolding transformation of each energy level {e n } — > {e n } is done by using the 
cumulative mean number of levels up to the energy e [20] , 




(Al) 
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In the above equation, d(e) exhibits the density of energy levels obtained by the Thomas- 
Fermi rule [5]: 

d(e) = V(e) =f6(e- H(q, p)) d'qd'p, (A2) 

where (2nK)f is the Planck volume of the system with / degrees-of-freedom, V(e) is the 
phase volume on the energy surface, 8(e — H) is the delta function, ff(q, p) is the classical 
Hamiltonian function, and (q, p) are the coordinates and momenta in the phase space. The 
unfolding transformation of spectrum {e Tt } — > {e n }(n = 1, 2, 3 • • •) is defined by 

e n = #([0,e n ]). (A3) 

Here, we consider the unfolding procedure of the i-th sub-spectrum. Since the phase 
volume of the i-th component is pi(e)V(e), the density of each sub-spectrum is then described 
by 

di(e) = = Pi (e)d(e) (A4) 

When the energy interval [e, e + Ae] is sufficient small, the phase space geometry on the 
energy surface does not change in general. In other words, Pi(e) is approximated by a 
constant value in this ivterval, and the cumulative mean number of the i-th sub-spectrum 
is thus described as 

Me,e + Ae]) = p t (e)V(e)de (A5) 

- (2^7 le +At V(e)de = Pl *^ e + AC]) - (A6) 



Therefore, in the asymptotic limit Ae — > 0, one can see that each spectral component obeys 
a same scaling law and is unfolded automatically by an overall unfolding transformation: 

en = 4{e n ) = -#i(en). (A7) 

pi 

The billiard system is very convenient since the phase space geometry does not change 
for variety of e, and #i(e) = Pi#(e) is satisfactory even when the size of energy interval does 
not small. 
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FIGURES 

(a) (b) (c) 
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FIG.l Numerical results of the level spacing distribution P(S) for various values of a ; 
(a) a = 1 + f x 1(T 4 , (b) a = 1 + f x 1(T 9 , and (c) a = 1 + ttx 10 -11 . We used energy levels 
e m>i e [300 x 10 7 , 301 x 10 7 ]. Total numbers of levels are (a) 10000016, (b) 10000046, and 
(c) 10000043. The dotted curve in each figure shows the Poisson distribution, P(S) = e~ s . 
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FIG.2 Function | log (1 - M(S; N)) \ for (a) a = 1 + f x 10~ 4 , (b) a = 1 + f x 10" 9 , and 
(c) a = 1 + 7r x 10~ n . The dotted line in each figure corresponds to the Poisson distribution, 
M(S) = l-exp(-S'). 
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FIG.3 D(S;N) for TV = 11339 and for N = 61905, which correspond to energy levels 
e m>i G [100 x 10 6 , 101 x 10 6 ] and e„ M G [300 x 10 7 ,301 x 10 7 ], respectively. In each case, we 
fixed a = 1 + | x 10~ 4 . The dashed line, D(S) = 1, exhibits the Poisson distribution. 
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FIG.4 D(S;N) for (a) a = 1 + f x 10~ 4 and N = 61905, (b) a = 1 + f x 10~ 9 and 
iV = 61906, and (c) a = l+7rx 10~ n and iV = 61906, which are calculated from energy levels 
e m> j G [300 x 10 7 ,301 x 10 7 ]. The dashed line, D(S)=1, exhibits the Poisson distribution. 
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